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By means of density functional theory (DFT) and the generalized gradient approximation (GGA) 
we present a structural, electronic and magnetic study of FePt, CoPt, FeAu and FePd based Llo 
ordered cuboctahedral nanoparticles, with total numbers of atoms, Nt ot = 13, 55, 147. After a con- 
jugate gradient relaxation, the nanoparticles retain their Llo symmetry, but the small displacements 
of the atomic positions tune the electronic and magnetic properties. The value of the total magnetic 
moment stabilizes as the size increases. We also show that the Magnetic Anisotropy Energy (MAE) 
depends on the size as well as the position of the Fe-atomic planes in the clusters. We address 
the influence on the MAE of the surface shape, finding a small in-plane MAE for (Fe,Co)24Pt3i 
nanoparticles. 



I. INTRODUCTION 

Current uses of nanometer designed magnetic devices 
range from biomedical applications [IH3], catalysis [I], 
energy harvesting [5] to data storage [BHS]- Because of 
these diverse applications in nanotechnology [TQl E] , the 
development of nanostructured magnetic materials has 
become a highly active field. Focusing in the field of mag- 
netic recording, a huge amount of experimental and the- 
oretical work has been carried out during the last decade 
to seek novel approaches to construct advanced materials 
for ultrahigh density magnetic storage, with the aim of 
increasing the state-of-the-art beyond 1 Tbit/in 2 ( [81IT2T- 
Il5p. Most approaches are focused on thin films or multi- 
layers [TBHT5] and recently on slabs Q15] . However, during 
the last decade has emerged the possibility to use clus- 
ters deposited on surfaces [20H26] to increase the record- 
ing density. These clusters or nanoparticles (NPs) have 
properties different from those of bulk alloys due to their 
reduced surface atomic coordination. In particular, bi- 
nary 3d-5d NPs formed by transition metals (TM) such 
as Fe or Co together with 5d noble metals like Au or 
Pt allow the possibility to tune the magnetic properties 
based on an in-depth knowledge of their geometrical [2"71 - 
[29] and magnetic behavior [30ll32j . 

One vital physical quantity in magnetic recording is 
the magnetic anisotropy energy (MAE) of the storage 
medium. The MAE determines the tendency of the mag- 
netization to align along some specific axis in solids and 
clusters. As we have pointed out, binary magnetic NPs 
based on (Fe,Co)Pt are good candidates for novel mag- 
netic recording media, especially those phases chemically 
Ll ordered, where the value of the MAE is of order 
of 7x 10 7 erg/cm 3 ( 33 ). The trend to higher record- 
ing densities requires continuous reduction in the grain 
size while retaining large values of the factor KV/kT to 
avoid loss of recorded information due to the onset of 
superparamagnetic behavior [M] [33] . Following the Neel 
relaxation law [36] . the only way to reduce the size of 
the NPs avoiding this trouble is through a higher val- 
ues of the MAE. To control this magnetic energy, mainly 



determined by the spin-orbit coupling (SOC) [37], it is 
necessary to investigate the structure, the electronic and 
the magnetic behavior of these systems. For this pur- 
pose, Gambardclla et al showed experimental and the- 
oretically that when Co adatoms were deposited onto 
a Pt(lll) surface, they had a MAE of 9 meV/at aris- 
ing from the strong SOC induced by the Pt substrate 
and for a unquenched orbital moments |20j . In addition, 
they increased the number of Co atoms on the metal 
surface forming NPs that ranged from 3 up to 40 atoms. 
The results showed that smaller NPs exhibited a higher 
MAE. These results opened a route to understand and 
fabricate high density magnetic recording materials us- 
ing deposited NPs on surfaces. There are several exper- 
imental [331 GH] and theoretical [2"5H321 |4"DP5] studies 
regarding isolated NPs aiming to obtain the best mor- 
phologies and magnetic behavior covering monometal- 
lic NPs HHH2I], binary alloys [H IH S3 US] and even 
capped NPs [33] |33] . Gruner et al have carried out a to- 
tal energy study of a wide range of structures of various 
shapes and sizes for (Fe,Co)Pt NPs 28 as well as for 
Fe(Pd,Ni) [35] .It was found that the most energetically 
favored structures obtained were those of ordered multi- 
ply twinned icosahedra and decahedra shapes. Gruner et 
al have also obtained locally the magnetic moment (MM) 
for Fe 2 56B 2 96, with B = Ni, Pd, Pt, Ir, Au, and as we 
will see in the present work, the tendency to augment the 
Fe MM in the vicinity of the cluster surface obtained by 
Gruner is in good agreement with our results. 

To obtain the MAE using the framework of DFT im- 
plies a huge computational resource since a fully rela- 
tivistic (FR) and a full potential (FP) treatment becomes 
necessary. A widely used approximation to overcome the 
all-electron (full potential) problem and to do quick and 
accurate calculations is substitute the core electrons by 
a pseudopotential (PP) [SI 07] . Most of the codes that 
use the PP approximation use the scalar- relativistic (SR) 
corrections (Darwin term and mass velocity) but they are 
not sufficient to address the behavior of magnetic systems 
because the MAE is mostly controlled by the SOC. Also, 
the magnetization density vector can vary from point to 
point in space presenting a spin non collinearity. To over- 



2 



come this barrier we have used a fully relativistic pseu- 
dopotential (FR-PP) approach implemented recently in 
the SIESTA code [SI GJ5]. 

The (Fe,Co)Pt Llo based alloys have large uniaxial 
anisotropy because of the layered structure (see Fig. [I]). 
The purpose of this paper is to investigate the struc- 
tural, electronic and magnetic properties of (Fe,Co)Pt 
and Fe(Au,Pd) Llo cuboctahedral nanostructured NPs 
having the total number of atoms N tot = 13, 55 and 147, 
and to calculate the MAE using the above mentioned FR- 
PP scheme. It is shown that the energy surface can be- 
come complex, showing features beyond the simple uni- 
axial anisotropy. This demonstrates the importance of 
investigating the dependence of the total energy on the 
orientation of the magnetization axis as we will see in 

Sec. Imol 

The paper is structured as follows. In section[ll]we de- 
scribe briefly the theoretical tools to perform all the cal- 
culations as well as the kind of NPs studied in the present 
work. The importance of the structural relaxations will 
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be explained in III A The local magnetic moments and 
the density of states are described in subsections |IIIB| 
and |III C| respectively. The MAE and its separate con- 
tributions is discussed in Sec. Ill D[ Finally, Sec. IV ; 
marizes the main results. 



sum- 



II. COMPUTATIONAL DETAILS 

We have undertaken calculations of electronic struc- 
ture and magnetic anisotropy energies (MAE) by means 
of DFT using a recent fully relativistic (FR) imple- 
mentation [54] in the GREEN [56J [57] code employ- 
ing the SIESTA [55] framework. We use fully separa- 
ble Kleinmann-Bylander |58j and norm-conserving pseu- 
dopotentials (PP) of the Troulliers-Martins [59] type to 
describe the core electrons. Our DFT based calcula- 
tions have been performed within the generalized gra- 
dient approximation (GGA) for the exchange correla- 
tion (XC) potential following the Perdew, Burke, and 
Ernzerhof (PBE) version [SU]. To address the descrip- 
tion of magnetic systems, pseudocore (pc) corrections 
were used to include in the XC terms not only the va- 
lence charge density but also the core charge as Louie et 
al [55] pointed out. In order to ease the convergence 
of three center integrals with the size of the real space 
grid, p c (r) is replaced by a pseudo-core charge density, 
p pc (r), which equals the real core charge density beyond 
a given radius, r pc , while close to the nuclei it becomes a 
smooth function. The radius r pc should be chosen small 
enough to ensure that the overlap region between the va- 
lence and the core charges is fully taken into account. 
Based on previous studies of the binary alloys [51], we 
have choosen for the radius that equals the core and va- 
lence charge the values of r pc (Fe,Co) = 0.6 Bohrs and 
r pc (Pt,Au,Pd) = 1.0 Bohrs, ensuring that the overlap re- 
gion between the valence and the core charge is fully take 
into account. As basis set, we have employed double-zeta 
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FIG. 1. (Color online) Initial Llo cuboctahedral structures 
employed in the simulations. Green spheres depict the mag- 
netic atoms while turquoise ones are the non magnetic species. 
On the NP at the extreme right, the yellow lines show two 
kinds of surfaces (square and triangle) . The three axis in the 
second NP (N to t=55) represent the usual cartesian frame, X, 
Y, Z being the angles 9 and <j> for each one of them (0°,0°), 
(90°, 0°) and (90° ,90°), respectively. The distance between 
planes are marked in red by c parameter and the lattice con- 
stant by a. 



polarized (DZP) strictly localized numerical atomic or- 
bitals (AO). The confinement energy, E c , defined as the 
energy cost to confine the wave function within a given 
radius was set to 100 meV. The Fermi-Dirac distribu- 
tion was used to obtain the ocupation numbers and the 
electronic temperature was set to 50 meV. 

In the SR pseudopotential approximation, the Kohn- 
Sham Hamiltonian 
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is diagonal in spin space and collinear spin is as- 
sumed [51]. In equation Q T is the kinetic energy, 

ylocal ig the fully local lon 

g-ranged potential commonly 
set to the I = radial component of the PP, V KB is the 
Kleinmann-Bylander (KB) operator [55], Vjj the Hartree 
term and Vxc is the exchange-correlation operator. Just 
two of those terms depend on the spin projections (say 
along the z axis), a (=f, !): the KB term and the fi- 
nal (exchange correlation) term. In the collinear case, 
there is a common quantization axis for the whole sys- 
tem, and the charge density has two independent projec- 
tions, p^(r) and p^(r), parallel and antiparallel, respec- 
tively. However, in the FR-PP aproximation, off-diagonal 
spin terms appear in the Hamiltonian causing a mixture 
of spin components because the spin quantization axis 
varies from point to point in space -i.e. non collinear 
case. Consequently it was necessary to use the scheme 
developed by Kiibler et al [53J that will give the mixed 
components for XC potential (see Ref. [53] for details). 

From an ah initio point of view the MAE is defined as 
the difference in the total energy between easy and hard 
magnetization axis. It is common to fix the spin quanti- 
zation axis as the z direction. However, when the FR-PP 
approximation is used and we need the total energy in 
several directions it necessary to proceed in a different 
way, specifically to generalize the magnetization direction 
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to an arbitrary axis, 5 U , characterized by polar angles 
and <j>. The procedure will give us a new set of matrix 
elements as a function of 9 and </> angles for the KB term, 
Vffi. For the total energy calculations required to deter- 
mine the MAE we obtain self-consistency by means of the 
Hamiltonian instead of using the density matrix. To this 
end, in each iteration the Hamiltonian is obtained after a 
Pulay mixing [64] of the input and output Hamiltonian, 
H m , H out respectively. The criterion for a self-consistent 
solution is the requirement that input and output values 
differ by less than 1 meV. For each different set of angles, 
($',<//), we restart the self-consistent scheme using as in- 
put Hamiltonian the one output for the previous angles 
(6,0): 

H f,in = Jjout _ yKB + yKB (2) 

where the primes denote the matrices calculated for an- 
gles (&,((>'). 

The unit cell for Llo metallic based alloys consists of 
two fee cells displaced along the diagonal of the cube. 
The presence of two different kinds of atoms generates a 
vertical distortion so that its structure is defined by two 
quantities, the in-plane lattice parameter, a, and the out- 
of-plane constant, c. Prior to relaxation the NPs were 
constructed from their bulk jet phase forming a perfect 
Ll ordered cuboctahedron (see Fig. [I]). We have re- 
stricted our study to the so-called magic cluster sizes N tot 
= (10n 3 +15n 2 +lln+3)/3, where n is the number of ge- 
ometrical closed shells, being the total number of atoms 
for each species Nj,/ =(5n 3 +6n 2 +4n)/3 for magnetic (M) 
species -Fe and Co-, and Natm =(5n 3 +9n 2 +7n+3)/3 for 
non-magnetic (NM) species -Pt, Au and Pd. The initial 
lattice parameters, a, as well as the c/a ratios were cho- 
sen as their bulk experimental values |65j : &p e pt = 3.86 
A and (c/a)_FePt = 0.98; &p e p c i = 3.89 A and (c/a)_F e pd 
= 0.938; &FeAu = 4.08 A and (c/a,)p eAu = 0.939; a CoPt 
= 3.81 A and (c/a) Co pt = 0.968. 



III. RESULTS 

We have carried out a systematic study of bimetal- 
lic nanoclusters, concentrating on the magnetic proper- 
ties (spin and MAE). The calculations were made on fully 
relaxed structures produced using a conjugate gradient 
method. In the following sections we present the results 
of the magnetic property calculations, but first we con- 
sider the structures themselves as a basis for the inter- 
pretation of the magnetic properties. 



A. Conjugate gradient relaxations 

To carry out relaxation of the NP structures, we have 
employed the conjugate gradient (CG) method, minimiz- 
ing the forces between atoms until they were less than 
0.03 eV/A. The optimizations have been done at a spin 



polarized SR level, and just to address the calculations 
of MAE, spin moments and density of states (DOS) a 
FR-PP scheme was included. 
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FIG. 2. (Color online) Lorentzian broadening of the 
bond distances between magnetic (djvf-M) and non- 
magnetic (djvM-iVA/) atoms, thick green and thin blue lines, 
respectively. Each one of the two columns (A,B), from bot- 
tom to top, depict the distances when the number of total 
atoms, N tot , increase from 13 to 147, (A.l) to (A. 3) for FePt 
and (B.l) to (B.3) for CoPt. The three vertical lines show 
a (black dashed), c (red solid) and the first nearest neigh- 
bours nn (black solid) experimental lattice values in their bulk 
phases. The values are provided in the text. 

In Fig. [2] we show the evolution of the values of the 
lattice parameters a, c and the first nearest neighbor dis- 
tances, nn, after a CG relaxation of the (Fe,Co)Pt NPs. 
Although the NPs experience only a small reconstruc- 
tion, this is enough to change the magnetic properties 
substantially as we will show in section |IIID| The dis- 
persion in the nn values for Fe species (green thick lines 
in the Fig. [2](A.1-A.4)) around the bulk value is ±0.4 A. 
This means that the Fe atoms have experienced a regular 
distribution around nn. Regarding Pt atoms (blue thin 
lines), this range is 0.2 A less than those of Fe atoms 
and also the average displacements are smaller. As a re- 
sult, the Pt atoms are concentrated closer to their bulk 
positions than Fe. With respect to the mean distance be- 
tween planes for each species, it is interesting to note that 
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Fe planes are closer after a reduction of the bulk value by 
-0.6 A. On the contrary the distances between Pt planes 
are larger, increasing by a value of 0.6 A. In general, we 
can say that the magnetic species has a higher dispersion 
around its bulk lattice parameters than the non magnetic 
one, except for the case of N tot = 55, where Pt atoms are 
also significantly distributed around c and a. It can be 
seen that for CoPt NPs (Fig. |j (B.1-B4)) the dispersion 
around a, c and nn is less than for FePt NPs. In this 
case, both Co and Pt atom positions deviate by ±0.15 A 
from their bulk nn structure values. As in the FePt case 
the Co atoms have reduced their mean separation values 
whilst those of the Pt atoms have increased. The distance 
between planes differ by smaller amounts than for FePt, 
the ranges being between -0.2 A and +0.3 A for M and 
NM atoms. This implies that for CoPt NPs there is less 
distortion of the bulk structure. The bond distances for 
Fe(Au,Pd) relaxed NPs (not shown here) have a similar 
behavior for each atomic species (magnetic and non mag- 
netic). Specifically, the Au atoms experience an increase 
in their nn distances of 0.15 A while the separation of Fe 
atoms decreases by 0.3 A. The out of plane variations are 
between +0.4 A for Au atoms and —0.3 A for Fe atoms. 
It is interesting to point out that in general the distance 
of the surface atoms from the center of the NPs tends to 
be reduced in comparison with the initial bulk structures. 
Studying this distances for the atoms located at differ- 
ent type of surfaces (squares on the top and the bottom 
and triangles or squares in the side of the molecules) (See 
Fig. [I]) , we can say that there is not a general trend either 
for magnetic nor non-magnetic species. 



B. Magnetic moments 

As a result of a Mullikcn analysis we show in Fig. [3] 
the variation of the spin magnetic moment (MM) values 
for every atom belonging to (Fe,Co)Pt NPs in function 
of its distance from the center. The MM values have 
been derived as the difference between the majority and 
minority spin charges and, although this is a qualitative 
study, it allows us to determine whether the NPs have 
more polarization in the surface or in the core. 

The local MM of the Fe and Co atoms (green filled 
squares), are remarkably large in comparison with those 
of the Pt atoms (blue empty squares). This is a general 
trend in all the clusters and the average differences range 
from 1.8 to 3.3 \ib for Co67Ptso and Fe67Ptso, respec- 
tively. Taking into account the region where the atoms 
reside (core or surface), the values will be slightly dif- 
ferent. So, we can observe that for magnetic atoms all 
the NPs have their inner MM values lower than those 
in the surface by ~0.4 hb- However, the Pt MM values 
remain around 0.25 /ib along their radial positions, the 
difference being only 0.1 \xb from inner to surface. This 
behavior prevails even for CoPt NPs. It is noticeable that 
the Co local MM are ~1.4 fis smaller than its magnetic 
counterpart in any other NP, even for Fe(Au,Pd) (not 
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FIG. 3. (Color online) Spin magnetic moments for each kind 
of atom (magnetic and non-magnetic) as a function of their 
distance from the center of the NP. Each row show, from left 
to right, how the MM values evolve when the total number 
of atoms increase for a specific alloy (Ntot = 13, 55, 147). Fe 
and Co atoms are depicted by green filled squares while Pt 
by empty blue ones. The red symbols show the MM values 
for each one of the atoms when Ntot = 13. Every graph has 
been divided by a solid line showing two zones that point out 
the core and surface regions. The added red dashed line in 
the first column also divide the smaller NP in the core and 
surface regions. 



shown in Fig. |3| . It is also interesting to note that the 
local MM at the surface in the smaller NPs (red symbols) 
have the largest values. The ratio of surface to volume 
atoms in these tiny NPs is 12/1 and 92% of the atoms are 
located on the surface. So, the surface effects are more 
pronounced at these sizes as we see in the increasing val- 
ues of the MM. 

In table U we summarize the total MM of all the NPs 
studied in this work. One of the main results is that 
for all the sizes, the FePt NPs have higher MM tot /N tot 
when compared with any other kind of NPs. The FePt 
values range from 1.62 /zs/at for FesPt 8 to 1.66 /is/at 
for Fe24Pt3i. Despite this small increase for MM tot /N tot , 
the different kind of NPs follow the same trend as the 
NP size increases. Fe67Pdso is an exception having a 
value of 0.15 /is /at less than Fe24Pd3i. If we inspect the 
third and the fourth columns, we note that the above 
mentioned increase of the MM jof /N tof is followed by the 
non-magnetic atoms but that the converse is true for the 
magnetic atoms. This loss of MM foi /Nji/ for the mag- 
netic atoms as the size of the NPs increases could be due 
to the fact that the percentage of surface atoms decreases 
from 80% with increasing N tot and as we have seen in 
Fig [3] the contribution of the higher spin values of the 
surface atoms will be diminished. The magnetic atoms 
are not entirely responsibility for the overall magnetic be- 
havior, the contribution from the spins of non-magnetic 
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atoms is vital to this complicated magnetic process. 

C. Density of states (DOS) 

To gain further insight about the electronic behavior 
of the NPs we present in Fig. [4] the spin resolved den- 
sity of states (DOS) projected onto Fe, Co and Pt atoms 
for both FePt and CoPt NPs, left and right panels, re- 
spectively. The atoms have been divided into two groups 
as previously in this work: surface (thick green lines and 
thin blue lines) and core (filled colored curves). The black 
lines show the total DOS. 

Firstly we note that for both types of NPs, as N tot in- 
creases, the total DOS peaks are smeared and bands are 
formed implying that the electrons become delocalized. 
Further, the projected DOS on magnetic surface atoms, 
together with their core counterpart, show that they pro- 
vide the largest contribution to the total MM of the NPs. 
The majority and minority bands of the smaller FePt and 
CoPt type NPs (A.l,B.l) have peaks around 0.25 eV and 
as the size of the NPs increases some of these peaks move 
below the Fermi level. For minority states, the FePt peak 
moves up to -0.3 eV, and up to the Fermi level for the 
majority states. CoPt NPs have the same behavior only 
for majority states while the minority peak remains at 
energies greater than for larger sizes. These displace- 
ments imply that the d bands are filling and as a result, 
there is a decrease of the total surface MM for the mag- 
netic species of 0.4 for FePt and 0.15 hb for CoPt (see 
first row in the Fig. [5]) as N tot increases. Other features 
of the total majority DOS of smaller NPs are the humps 
at -1 eV, -4 eV and at -5 eV for FePt and at -1 eV, - 
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TABLE I. Magnetic moment (MM) values in /z fl /at of all Ll 
cuboctahedral based NPs. The first column displays the four 
kinds of alloyed NPs and the second shows the number of total 
atoms, N tot . The MM values in the third and fourth column 
have been calculated by means the equations MM(„( /Nm and 
MMtoi/NjvA/, respectively, where as in the text, M refers to 
magnetic atoms and NM to the non-magnetic ones. The last 
column show MM tot /Nt ot . 
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FIG. 4. (Color online) Spin resolved density of states (DOS) 
of FePt (A. 1-3) and CoPt (B.l-3) alloyed NPs. Full black lines 
show the total up and down DOS when the total number of 
atoms, Nt ot , varies from 13 to 147 from bottom to top. The 
projected DOS of the surface atoms are represented by thick 
solid green lines for Fe and Co atoms and thin solid blue lines 
for the Pt atoms. The filled curves show the projected DOS 
for the core M and NM atoms, green and blue, respectively. 



3.5 eV and -5 eV for CoPt. The first two peaks located 
at higher energies persist for larger NPs, however the last 
disappears when N tot >13. It is wort noting that the 
FePt NPs with N tot >13 have a dip in the minority DOS 
at the Fermi level, showing that the minority channel is 
dominated by Fe surface atoms. This is in good agree- 
ment with the work of Gruner et al. so that this feature 
gives a way to distinguish between different morpholo- 
gies such as icosahedron or Llo cuboctahedron [27]. The 
DOS of FePd and FeAu NPs (not shown here) present a 
slightly different shape, but nonetheless exhibit the main 
feature the primary responsibility of the M atoms for the 
polarization in these NPs. 



D. Magnetic anisotropy 

We finally present the calculations of the magnetic 
anisotropy (MAE). In order to get a better knowledge 
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of the magnetic behaviour of the NPs, we also show in 
Fig. [5] together with the MAE, the MM for Fe, Co, Pt, 
Pd and Au atoms in the first and second rows. It is easy 
to distinguish between the MM of the surface and core 
atoms, whether or not they are magnetic, since as we 
have seen in sec. IIIIBI that the local surface MM values 
are higher than those of the core. 

Although there are some common tendencies in the be- 
havior, there is no overall trend, presumably because of 
the complexity of the atomic rearrangements and charge 
transfer. Consider first the behavior of the MM values 
of the magnetic atoms in FePt, CoPt, FePd and FeAu. 
The common factor in the behavior of all systems is an 
increase of the MM in the surface over that of the core 
atoms. In addition, FePt, CoPt and FePd exhibit large 
differences (as large as Afi = 0.7 /Us/at for FePt), which 
decreases with increasing N tot . The similarity presum- 
ably reflects the chemical similarity of Pd and Pt. Al- 
though its surface atoms have a larger MM than the core 
atoms, FeAu breaks the trend in that A/x remains reason- 
ably constant, presumably reflecting the different atomic 
rearrangements and charge transfer. Turning to the MM 
values of the non magnetic species, the tendency of the 
MM is to be almost constant within both the core and 
surface regions. Again, we note that the non magnetic 
atoms of FeAu NPs exhibit a different trend, and that 
further their MM values for N tot = 55, 147 are negative. 
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The calculated values of the MAE are shown in 
Fig. [5]JA.l,B.l). The alloys from which our NPs have 
been constructed have in their bulk Llo phases a MAE 
of order of a few meV [6 5) with the easy magnetization 
axis perpendicular to the atomic planes (see Fig(T]). We 
will see that most of all the studied NPs have the same 
easy axis orientation as their bulk alloys. Also, the val- 
ues of the total MAE obtained in this work are of order 
of tens of meV following the same trend shown by other 
groups for small NPs [3H [33- In the graphs, the MAE 
is expressed in meV per atom by dividing by the total 
number of atoms (magnetic plus non-magnetic) of each 
NP and using straight coloured lines we show the MAE 
values for each Ll alloy. 

Consider first the case of FePt and CoPt shown in 
Fig. |] (A. 1). Although we do not have site-resolved MAE 
values, we can interpret the data in relation to previous 
calculations of bulk properties of FePt [tJrJl RTF] . These 
suggest that the primary contribution to the MAE in 
FePt is a 2-ion anisotropy of the Fe sites mediated by the 
Pt sites. This suggests that the presence of surfaces and 
the consequent loss of coordination might be expected to 
lower the overall MAE, which is certainly the case for 
the two larger NP sizes considered here. However, it is 
interesting to note that the smallest NP size exhibits an 
increased MAE. Although we cannot here obtain site- 
resolved information for the MAE, it seems reasonable 
to suggest that this arises from the modified electronic 
properties within the smallest NPs. This is worth fur- 
ther consideration, with site-resolved calculations, since 



FIG. 5. (Color online) Magnetic anisotropy (MAE) values per 
atom and mean surface and core magnetic moments (MM) 
values per atom for Fe, Co, Au, Pd and Pt of FePt (green), 
CoPt (blue), FeAu (red) and FePd (turquoise) NPs as funtion 
of their total number of atoms. 



this enhanced MAE may be useful for applications. 

Regarding the orientation of the easy axis, most of the 
NPs studied present easy axis along the Z-axis. However, 
we note that for Fe 5 Pd 8 , Fe 5 Au 8 and Fe 6 7Pd 8 o the MAE 
has a negative value which means that the easy axis lies in 
the XY-plane. Further evidence of contributions to the 
MAE beyond the simple uniaxial case is shown in Fig. [6] 
Here we show the variation of the total energy of the 
Fe24Pt3i (upper row) and Co24Pt3i (lower row) NPs with 
the 9 (left) and <j> (right) angles. In both types of NPs 
the easy magnetization axis lies along the (001) direction 
-having the minimum value of the energy (9 = 0° , <f> = 
0°). Fixing (f> to 0° (empty blue squares) and 45° (full 
green dots) and varying 9 from zero to 180° we obtained 
different maxima for Co24Pt3i while the Fe24Pt3i NPs 
exhibit purely uniaxial behavior, with no dependence of 
4>. The graphs on the right side sweep the energy from 
4> = 0° to <fi = 180° keeping 9 constant. It can be seen 
that the in-plane magnetization for CoPt has two minima 
exactly at 45° and at 135° degrees (see Fig. [I]). In the 
case of FePt NPs, no in-plane anisotropy is observed. 
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FIG. 6. (Color online) Energy variation as a function of the 
f5(left) and </>(right) angles for Fe24Pt3i and Co24Pt3i NPs in 
the upper and lower rows, respectively. The zero of energy is 
set to the minimum value of E and all the points have been 
joined with lines in order to guide the eye. 



IV. CONCLUSIONS 

In conclusion, we have carried out a fully relativistic 
calculations, within the GGA approximation, of the mag- 
netic moments, density of states and MAE of Llo cuboc- 
tahedral FePt, CoPt, FeAu and FePd based NPs. We 
have restricted the total number of atoms to the magic 
numbers: 13, 55 and 147, giving diameters of the NPs 
from 0.6 nm for N tot =13 up to 1.6 nm for N tot =147. 

Although the original stacking is retained after CG 
relaxation, the atoms exhibit small displacements from 
their original bulk positions. The bond distances be- 
tween magnetic species have higher dispersion around 
the characteristic bulk values than exhibited by non mag- 
netic atoms, the non magnetic species being almost at the 
same position. Although this trend is followed by most 



of the NPs, there is an exception for CoPt NPs that show 
much less dispersion for both magnetic and non magnetic 
species. 

Regarding the magnetic structure we have shown that 
the outermost local magnetic moments of all the NPs 
studied in this work are larger than in the core, in good 
agreement with previous investigations. This magnetic 
behaviour is correlated with the PDOS analysis that 
shows the importance of the magnetic ordering of the 
surface magnetic species polarization. Also we showed 
that the MAE is size and stacking dependent and that 
the value increases for the smallest NPs respect to the 
bulk values. This indicates enhanced thermal stability of 
the smallest NPs. However, the larger FePt and CoPt 
NPs showed a reduction of the MAE consistent with the 
loss of coordination at the surface and a consequent re- 
duction of the (dominant) 2-ion anisotropy. This is an 
interesting observation which shows a dramatic change in 
the magnetic behavior in the smallest NPs which is worth 
investigating using site-resolved MAE calculations. 

The easy magnetization axis generally lies along the 
(001) direction, although in some FeAu and FePd NPs 
the anisotropy lies in-plane. As an example of an in- 
plane magnetic anisotropy we obtained 9 and 4> energy 
dependence for (Fe,Co)24Pt3i NPs in the Fig. [6] showing 
that the surface shape it is important to study the entire 
energy surface to investigate the overall form of the MAE, 
which, certainly for the case of Co24Pt3i has a significant 
contribution from a cubic anisotropy term in addition to 
the main uniaxial term. 
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